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Topological phases like topological insulators or superconductors are fascinating quantum states of matter, 
featuring novel properties such as emergent chiral edge states or Majorana fermions with non-Abelian braiding 
statistics. The recent experimental implementation of optical lattices with highly tunable geometry in cold gases 
opens up a new thrust on exploring these novel quantum states. Here we report that the topological non-trivial 
Bloch bands can arise naturally in a noncentrosymmetric lattice. It induces a controllable orbital hybridization, 
producing the topological band structure. In such bands, when considering attractive fermionic atoms, we 
find a topological Fulde-Ferrell superfiuid state with finite center-of-mass momentum in the presence of onsite 
rotation. This topological superfiuid supports Majorana fermions on its edges. Experimental signatures are 
predicted for cold gases in radio-frequency spectroscopy. 


Symmetry plays an important role in solid state materials 
and influences many of their properties in a profound way. Re¬ 
cently, the noncentrosymmetric materials, i.e., crystal struc¬ 
ture lacking a center of inversion, have attracted considerable 
theoretical and experimental attention in various fields of con¬ 
densed matter physics |[ll] . In contrast to the centrosymmetric 
case, the absence of inversion symmetry leads to very rich 
physical phenomena, such as skyrmion states IM , novel su¬ 
perconducting phases |[1-[1], as well as magnetoelectric ef- 
feet d m. In parallel to the developments in solid state 
physics, optical lattices with highl y tu nable geometry in the 
recent ultracold atom experiments Bill] provide new opportu¬ 
nities to study noncentrosymmetric materials. Taking the ad¬ 
vantage of high controllability, the cold gases are expected to 
emerge not only as a flexible tool to simulate electronic sys¬ 
tem, but also a versatile platform to create fascinating quan¬ 
tum states of matter not existing in nature. 

Recent experimental progress in manipulating higher or¬ 
bital bands in optical lattices ifT^Ul^ provides unprecedented 
opportunities to investigate quantum many-body phases with 
orbital degrees of freedom. Studying higher orbital physics 
in optical lattices is attracting considerable interests due to 
their unique and intrinsic spatial nature d. An idea of using 
orbital hybridization to emulate spin-orbit physics, or artifi¬ 
cial gauge fields in general, has emerged in recent theoretical 
studies iflSl - l^ . where various interesting many-body states 
have been proposed. In the context of periodically driven sys¬ 
tems, orbital hybridization is shown to be controllable by lat¬ 
tice shaking ll23l - [^ . However how to systematically control 
the orbital hybridization in static optical lattices remains un¬ 
clear and stands as an obstacle to explore the rich phenomena 
in orbital hybridized many-body ground states. 

In this letter, we propose a static noncentrosymmetric opti¬ 
cal lattice, where the orbital hybridization is systematically 
controllable by manipulating symmetry breaking. We find 
topological bands arise from the interplay between high or¬ 
bitals and inversion symmetry breaking, yet without requir¬ 
ing Raman induced spin-orbit coupling nor other artificial 
gauge field inm. Furthermore when considering attrac¬ 
tive fermionic atoms, say ^Li ll28l - [^ loaded into such a non¬ 


centrosymmetric optical lattice, we find an orbital hybridized 
topological Fulde-Ferrell superfiuid state (tFF) in the pres¬ 
ence of local orbital angular momentum induced by onsite 
rotation |[^ . This state represents an unconventional type 
of Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state, topolog¬ 
ically distinct from the 5-wave FFLO superfluidity observed 
in a quasi-one dimensional Fermi gas ll^ . This orbital re¬ 
alization of the tFF superfiuid state here is distinguished from 
previous works based on hyperfine states ilssi - l^ . where spin- 
orbit coupling and Zeeman fields are required. Moreover, 
our symmetry-based method of controlling orbital hybridiza¬ 
tion is in principle applicable to other optical lattice setups 
as well |[l8l - [^[^ l37h. which could lead to more interesting 
noncentrosymmetric many-body phases worth future explor¬ 
ing. 

Effective model Let us consider a noncentrosymmetric op¬ 
tical lattice with the potential 

V{x,y) = -Vx cos^(kLxx) - Vy cos^ (kLyU) 

+ VYCosy3kLyy + e/2), ( 1 ) 

where Vx, Vy, Vy are the depth of optical lattices, kpx, kpy 
are the wavevectors of laser fields and the corresponding lat¬ 
tice constants are = njkLx, cty = Trlkpy along x and y 
directions respectively. By the techniques of designing the ge¬ 
ometry of optical lattices developed in the recent experimental 
advances UU\ , the configuration of optical lattices considered 
here can be formed through three retro-reflected laser beams 
as shown in Fig. (Ha). The interference of two perpendicular 
beams X and Y gives rise to a 2D square lattice. The third 
beam Y creates an additional standing wave pattern which 
breaks the inversion symmetry along the y direction, for ex¬ 
ample, when 0 = shown in Fig. [TJb). This noncen¬ 

trosymmetric geometry plays a crucial role in producing the 
non-trivial Bloch bands in our model, to be illustrated below. 
Here we consider highly anisotropic geometry of the optical 
lattice through increasing lattice depth (Vy ^ Vx) and spac¬ 
ing (Oy ^ Ox) in the y direction, meanwhile requiring that 
the lattice potential satisfy the condition Vxkj^^ = Vykj^y. 
This should keep local rotational symmetry of each site in the 
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FIG. 1: Noncentrosymmetric optical lattice and topological band 
structures, (a) Three retro-reflected laser beams create the lattice po¬ 
tential in Eq. Q. X and Y interfere and produce a 2D square lattice, 
while Y creates an independent standing wave, (b) Lattice potential 
along the y direction shows the inversion symmetry breaking, (c) 
The single-particle energy spectrum along kx axis in the unit of Erx 
for the lowest four bands through plane wave expansion calculation. 
Here the lattice depth is chosen as Vx / Erx = 4, Vr / ERy = 40 
and Vy!ERy — 4 with the recoil energy Erx = h‘^k‘l^/2m and 
ERy = h^kly/2m. 

xy plane ll^ . The band structure of such a lattice system is 
solved numerically through plane wave expansion. We find 
that the second and third bands cross in the absence of Vy 
with the band touching points protected by time-reversal and 
inversion symmetries. With finite Vy, we see the gap reopen¬ 
ing due to inversion symmetry breaking (Fig.[TJc)). The band 
mixing can thus be turned on and off by controlling the sym¬ 
metry of the lattice geometry. 

The essential physics of band mixing is captured by the fol¬ 
lowing multi-orbital tight binding model without considering 
the full band structure theory, 

Ho = t^^Cl^{Y)Cp^{Y + e:c)-tyY(^lyi^)Cpyi^ + ^^) 

r r 

+ (r + 4) - i XI Cl^ (r)Cp^ (r + 4) 

r r 

- + 4X 

r r 

+ /i.c.+ C^^(r)C'p^(r)], (2) 

r 

where Cjy (r) is a fermionic annihilation operator for the lo¬ 
calized u orbital (px or py) located at the lattice site r and 
the chemical potential is denoted by fi. Since the system is 
relatively stronger confined in the y direction, the tunnelings 
in the y direction are much weaker compared to that 
(tx^ ty) along the x direction. The relative sign of the hopping 


FIG. 2: Bulk and edge properties of the topological band insulator. 

(a) The single-particle band structure calculated from the tight bind¬ 
ing model (Eq. (|2])) in quasi-lD limit, when tx/ty = 8, t/ty = 0.6. 

(b) Topological winding of Hamiltonian in Eq. © across the Bril- 
louin zone, (c) Energy spectrum of the Hamiltonian (Eq. (IJJl) in 
quasi-ID limit with finite length (lattice sites N — 200), when 
tx/ty = 8 and t/ty = 0.9. There are two zero-energy states in¬ 
side the gap which are located at the two outer edges of the system 
respectively. Here n labels the energy level, (d) The effect of the 
small transverse tunneling. The midgap bands show a finite curvature 
along the y direction with open (and periodic) boundary conditions 
in the x (y) directions when tx/ty = S, t/ty = 0.1, ty/ty = 0.02 
and t'x/ty = 0.1. 

amplitude is fixed by the parity of Px and Py orbitals. The key 
ingredient in our model is the hybridization between px and Py 
orbitals. It arises from the asymmetric shape of the Py orbital 
wavefunction induced by the inversion symmetry breaking in 
the y direction. This asymmetry leads to the orbital hybridiza- 
liontY,r[ClJr)Cp^{Y + e^) - Ct^(r)Cp,(r + e^)] + h.c. 
in Eq. Q), which plays a crucial role in producing topolog¬ 
ical non-trivial band structures. It resembles spin-orbit cou¬ 
pling |[38| - [4Q|] when the Px and Py orbitals are mapped to 
pseudo-spin-1/2 states. But unlike the spin-orbit coupling, 
the orbital hybridization is geometrically controllable through 
manipulating the inversion symmetry of the optical lattice. 

Topological band structures. We first focus on the quasi- 
one dimensional limit where the transverse {y direction) tun¬ 
nelings are negligible. At half filling (one particle per lattice 
site), in the basis of (k), (k)), the Hamiltonian takes 

a suggestive form 

i7o(k) = /io(k)I + h(k).(7, (3) 

where ho{k) = {tx - ty) cos{kxax), hx(k) = 0, hy{k) = 
-2tsm{kxax), hz{k) = {tx + ty) cos{kxax), and cFx,y,z are 

Pauli matrices. The energy spectrum reads E± = /io(k) =b 

h‘^{k) + h‘l{k). As shown in Fig. [2a), the system is in an 
insulating state with a band gap determined by the orbital hy- 
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bridization t. It is a Z 2 topological insulator protected by time 
reversal and reflection symmetries El . To visualize the topo¬ 
logical properties of the band structures, we show the vector 
h(k) winds an angle of 27r when the momentum k varies from 
—TT to TT crossing the entire Brillouin zone (BZ) in Fig. [2b). 
It is also confirmed from the calculation of winding number 
defined as W = f with h = -^ and 

^yz = —^zy = 1 - When t 7 ^ 0, the winding number is 1, 
which signifies a topological band insulator state. The non¬ 
trivial topology of this state also manifests through the exis¬ 
tence of the edge states. As shown in Fig. [2c), there are two 
emergent zero-energy modes located at the two outer edges of 
the system respectively. 

Next we will discuss the effect of the small transverse (y di¬ 
rection) tunneling, which has been neglected above but always 
exists in a realistic quasi-one dimensional system. By consid¬ 
ering small transverse tunneling, the zero-energy modes of in¬ 
dividual chain will morph into a midgap band, with finite cur¬ 
vature in the transverse direction as shown in Fig. [2d). The 
topological band insulator state remains stable at small value 
of transverse tunneling. For example, when tx/ty = 8 and 
t/ty = 0.1, the topological state survives until ty/ty reaches 
0.073 with = O.lty. However, beyond this value, the 
band gap will close and the topological band insulator state 
becomes unstable. 

Topological Fulde-Ferrell state and Majorana fermions. In 
this section, we study attractive fermions in the topological 
noncentrosymmetric optical lattice (Eq. ([T])) with onsite rota¬ 
tion |[^ and show that a topological Fulde-Ferrell superfluid 
state (tFF) with finite center-of-mass momentum emerges. 
The interacting model to describe this fermionic system is 

H = + Hint + Hl + Hz‘ (4) 

In this model, the interaction Hint = 

UY.rCl^\Y)Cp^{v)Cl^{Y)Cp^{Y) can be indeed by 
considering optical Feshbach resonance [^, Bose- 
Fermi mixtures El or dipolar atoms\molecules El. 
Hl = i^z WC'px (r)] is the onsite 

rotation induced orbital Zeeman energy. Hz the orbital 
splitting hZ,r[CUr)Cp^{Y) - CI^{y)Cp^{y)] originated 
from the onsite energy difference between px and Py orbitals, 
which is tunable through adjusting the lattice depth or lattice 
spacing in x and y directions. In the presence of onsite 
rotation, i.e., Qz 7 ^ 0, the Fermi surface of the system 
becomes asymmetric along the x direction. The pairing with 
the center-of-mass momentum Q and — Q are no longer 
degenerate. The Fulde-Ferrell (FF) state with a plane wave 
order parameter oc is more energetically favorable as 

compared to the Larkin-Ovchinnikov state with a stripe order 
oc cos(Q • r). 

Taking the superfluid pairing order parameter A(r) = 
U{Cpy (r)Cp^ (r)) = the system is described by the 

Bogoliubov-de-Genes (BdG) Hamiltonian at mean held level. 
Since the weak transverse hopping introduces a small Fermi 
surface curvature, we expect the center-of-mass momentum 
of pairing Q pointing along the x direction, say Q = Q(l, 0), 
in order to maximize the phase space of pairing. Through 


diagonalizing the BdG Hamiltonian, we obtain the free en¬ 
ergy of the system. The pairing order parameter A and the 
center-of-mass momentum of pairing Q are determined from 
minimizing the free energy. The details are given in the Sup¬ 
plemental Material. We And that when (2;^ > 0 the center-of- 
mass momentum of pairing along the x direction is selected 
as a certain positive Q due to the deformation of the Fermi 
surface. Since the Px and Py bands have different bandwidths, 
we introduce two dimensionless quantities flp^ = and 
jipy = which respectively control the average Ailing of 
the two p bands. Correspondingly, we also introduce a dimen¬ 
sionless chemical potential jl = and orbital polariza¬ 
tion h = Por Axed jl, the resulting phase diagram 

as a function of h and the orbital hybridization t is shown in 
Fig. [2a). There are two first order phase transitions as the po¬ 
larization h is increased (except dLitjty ^ 0.8). The first one 
is a transition from a gapped FF superfluid to a gapless FF 
superfluid state (with Bogoliubov quasi-particles being gap¬ 
less). The gap closing across the phase transition is shown 
in Fig. [2b). Further increasing the polarization, the second 
one occurs between the gapless FF and tFF superfluid states. 
When t/ty ^ 0.8, there is only one phase transition from FF 
to tFF superfluid states without passing through the gapless 
FF superfluids, since we And in Fig. [2c) that Eg will firstly 
close and reopen immediately. We also And that a finite po¬ 
larization is required to stabilize the tFF superfluid state. The 
critical polarization he decreases as the orbital hybridization t 
increases. 

The transition from non-topological to topological FF states 
here can be understood by observing the non-trivial Z 2 topo¬ 
logical invariant. The BdG Hamiltonian (See details in the 
Supplemental Material) maintains the particle-hole symme- 

try, i.e., SiJBciG(k)S-i = with ^ 

while the time reversal and chiral symmetries are broken. 
Therefore, the tFF superfluid state predicated here belongs 
to the D symmetry class according to the general classifica¬ 
tion scheme of topological superconductors j^ . This topo¬ 
logical state is thus characterized by a Z 2 topological invari¬ 
ant in, 113]. As shown in Fig. [2b) and (c), we And that the 
Z 2 topological invariant M = — 1 in the tFF superfluid state 
(see details in the Supplemental Material). 

To further demonstrate the topological nature of the tFF 
superfluid phase, we will show Majorana fermions are sup¬ 
ported in this state. To see this, we consider a cylinder geom¬ 
etry of the system, where the open (periodic) boundary con¬ 
dition is chosen in the x{y) direction respectively. The en¬ 
ergy spectrum in Fig. [2a) is labeled by the momentum ky. 
As shown in Fig. [2a), all the bulk modes are gapped and 
there are two degenerate flat bands composed of Majorana 
fermions located at the two outer edges of the system respec¬ 
tively. As shown in Fig. [2c), for a fixed ky, there are two 
zero-energy states located at the two outer edges of the sys¬ 
tem. The corresponding wavefunctions satisfy 

the relation = “'^0 c)n the 

right [left] edge (Fig. 2d) and (e)). These eigenstates support 
localized Majorana fermions at the edges of the system. These 
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FIG. 3: Zero-temperature phase diagram with fixed jl. (a) When the 
polarization h is small, the FF superfluid state is the ground state 
of the system. By increasing h, a first order phase transition occurs 
between FF superfluidity and a gapless FF superfluid state. Further 
increasing h, the system eventually evolves into topological FF su¬ 
perfluid state (tFF). Other parameters are tx/ty = 8, ty/ty = 0.05, 
t'x/ty = 0.1, U/ty = —15, jl — 0.4, and VLzjty — 0.3. (b) and 
(c) The quasi-particle exaction gap Eg as a function of polarization 
h when tjty = 0.6 and tjty = 0.8 respectively. The vertical lines 
mark the point where the Z 2 topological invariant changes. 


Majorana fermions are signified through the local density of 
states (LDOS) which can be measured by radio-frequency (rf) 
spectroscopy lld^l - l^ . The LDOS is calculated sls p{x,E) = 
1/2 E^v{Wn,u\‘^S{E-(n) + \vn,u\‘^S{E + (ri)], where 
{un,iy, Vn,vY' is the eigenvector corresponding to the eigenen- 
ergy (jn of the mean-field BdG Hamiltonian with cylinder ge¬ 
ometry. We find that the zero-energy Majorana fermions man¬ 
ifest themselves by a peak in LDOS located at the edges of 
the system, as shown in Fig. Ob). This spatially localized 
zero-energy peak in LDOS can be detected using spatially 
resolved radio-frequency (rf) spectroscopy technique iiai, 
which would provide a concrete signature for the experiment. 

We also consider the effect of small transverse tunneling on 
the stability of tFF superfluids. The transverse hopping term 
restores the 2D nature of the Fermi surface, which acquires a 
flnite curvature in the transverse direction. It suppresses the 
perfect nesting, therefore it disfavors the tFF superfluid state. 
Our numerical result indicates that the tFF superfluity remains 
stable at small value of transverse tunneling. For example, the 
tFF superfluid state, as shown in Fig.^a) with t/ty = 0.8 and 
h = 0.7, survives until t'y/ty reaches 0.2 and the flat bands 
composed of Majorana fermions in a cylinder geometry will 
be maintained 11^ . Beyond this value, the bulk gap will close 
and the tFF superfluids become unstable. 

We would like to emphasize that the topological non-trivial 
properties in the noncentrosymmetric optical lattice arise di- 



FIG. 4: Topologically protected edge states eomposed of Majo¬ 
rana fermions in the tFF superfluid state, (a) Energy spectrum 
of the Hamiltonian in Eq. under the mean-field approximation 
with open (and periodic) boundary conditions in the x (y) direc¬ 
tions. There is a doubly degenerate flat band eomposed of Majorana 
fermions when the transverse tunneling is small, (b) The loeal den¬ 
sity of states (LDOS) (see the main text). The zero-energy peaks of 
the LDOS are located at the two edges of the system respeetively. 
(c) Energy spectra with kytty = 0. There are two degenerate zero 
energy states, (d) and (e). The wavefunctions of the two Majorana 
zero-energy modes as shown in (c), which satisfy 
the right edge and = ~'^o,px (other compo¬ 

nents of these two wavefunctions, say and , also satisfy 
the same relations). These two states support two local Majorana 
fermions at two outer edges of the system respectively. Here we 
ehoose finite length along the x direction with lattice sites N = 400 
and h = 0.7, tjty — 0.6. Other parameters are the same as in Eig.[3] 
and n denotes the energy level. 


rectly from symmetry-breaking induced orbital hybridization. 
This plays a curial role in producing the topological FF su¬ 
perfluid state. We expect noncentrosymmetric optical lattices 
would provide a fertile ground to support orbital hybridized 
topological phases in general. 

Conclusion. We have established a systematic approach 
to control non-trivial orbital hybridization in a static noncen¬ 
trosymmetric optical lattice, which is shown to give rise to 
unconventional topological properties. We And that this lat¬ 
tice system in the presence of onsite rotation supports a novel 
topological FF superfluid state when loaded with attractive 
fermions. It features Majorana zero energy modes located at 
the outer edges, leading to signatures in the local density of 
states as a concrete experimental evidence of the topological 
superfluid state. 
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Supplementary Materials 

S-l. BDG EQUATION IN MOMENTUM SPACE 


Through introducing the superfluid pairing order parameter A{t) = U{Cp^ (r)Cp^ (r)) = the system can be described 

by the Bogoliubov-de-Genes (BdG) Hamiltonian at mean held level 

/^Px(Q/2 + k) e{Q/2 + k) 

HBdoO^) = \ 

V A* 0 

where the Nambu basis is chosen to be (Q/2 + k), (Q/2 + k), (Q/2 — k), (Q/2 — k))^, ^p^(k) = 

2 tx cos{kxax) — cos{kyay) — {ja — h), ^p^(k) = —2ty cos{kxax)^2ty cos{kyay) — and5(k) = 2 itsm{kxax)-\-i^z- 

Then, the free energy can be obtained by diagonalizing the BdG Hamiltonian in Eq. dST]) by standard procedure as 
F[A] = l/2^[e,^(k)+£,^(k) + ^e(-E,(k))i;,(k)]-^M, 

k A 

where Ex is the quasi-particle energy, and 0 is the Heaviside step function. The pairing order parameter A and the center-of- 
mass momentum of pairing Q can be determined from minimizing the free energy. 


^Pcc(Q/2 k) 


A 
0 

e{Q/2 - k) 


(SI) 


-e(Q/2 - k) -5. (Q/2 - k)y 


S-2. Z2 TOPOLOGICAL INVARIANT 

To characterize the topological nature of the tFF superfluid state, we calculate the Z 2 topological invariant. Here, we introduce 
the Majorana operators as 7 ^(r) = (r) + Cp^ (r), 7 ^(r) = [Cp^ (r) - (r)]/i, 7 '^(r) = (r) + Cp^ (r) and 7 ^(r) = 

[Cpy (r) — (r)]/*, which fulfill the relations 7 ^'^(r) = 7 “(r) and the anticommutation relations { 7 “, 7 ^} = 2 Sa/ 3 S{r — r') 

with a or /3 taking A, B, C or D. In terms of Majorana operators, the Hamiltonian in Eq. (4) under the mean-flied approximation 
can be represented as Hmf = Aiminm with 74^-_3 = 74 j -2 = 74 i-i = and 74 ^ = 

7 ^(rj), where j runs over all the N lattice sites and A is a skew-symmetric matrix. The Pfafflan of matrix A is deflned as 
Pf(A) = 2 ^^ SrG 52 nm=i ( 2 m-i),r( 2 m) with Ti = 2N, where 8271 IS the set of permutation and sgn{r) is the 

corresponding sign of that. The Z 2 topological invariant is deflned as M = sgn[Pf{A)] when choosing the periodic boundary 
condition. The Z 2 topological non-trivial phase is characterized by M = — 1, where as the topological trivial phase corresponds 
to M = 1. ._ 
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